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ABSTRACT 

We present the Lensed Mock Map Facility (LeMoMaF), a tool designed to perform 
mock weak lensing measurements on numerically simulated chunks of the universe. 
Coupling N-body simulations to a semi-analytical model of galaxy formation, LeMo- 
MaF can create realistic lensed images and mock catalogues of galaxies, at wavelengths 
ranging from the UV to the submm. To demonstrate the power of such a tool we com- 
pute predictions of the source-lens clustering effect on the convergence statistics, and 
quantify the impact of weak lensing on galaxy counts in two different filters. We find 
that the source-lens clustering effect skews the probability density function of the con- 
vergence towards low values, with an intensity which strongly depends on the redshift 
distribution of galaxies. On the other hand, the degree of enhancement or depletion in 
galaxy counts due to weak lensing is independent of the source-lens clustering effect. 
We discuss the impact on the two-points shear statistics to be measured by future mis- 
sions like SNAP and LSST. The source- lens clustering effect would bias the estimation 
of (78 from two point statistics by 2% — 5%. We conclude that accurate photometric 
redshifts for individual galaxies are necessary in order to quantify and isolate the 
source-lens clustering effect. 

Key words: cosmology: gravitational lensing - large-scale structure - methods: nu- 
merical 



1 INTRODUCTION 

One of the key goals of modern astrophysical and astro- 
nomical research is to map out the spatial distribution of 
the various matter components of the Universe. In the ac- 
tual concordance cosmological model, about 85% of the mat- 
ter in the universe is thought to be non-baryonic and non- 
interacting dark matter, the other 15% being composed of 
baryons (Spergel et al. 2006:). 

One of the most promising tools to track the dis- 
tribution of dark matter on cosmological scales is weak 
gravitational lensing, which has already proven to be well 
suited for other purposes, such as precision cosrnology 



suitea tor otner purposes, sucn as precision cosmology 
||HuI|2002|: iHuterer 20 02: Benabed fc Van Waerbekd 12004 : 
iBernstein fc Jainll20o4 Ixereno et al.ll2005l ). 

Weak gravitational lensing affects the observed galaxy 
properties such as ellipticities, magnitudes and apparent po- 
sitions in the sky. In the weak lensing regime, these effects 
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can only be measured in a statistical sense. Indeed, the 
first detections of the weak lensing sig n al (Bacon et al 200^; 
Kaiser et al] l2000l: iMaoli et all I2001I: IVan Waerbeke et all 



2000l . l2002l : lBrown et al.ll2003l : |jarvis et aLlbOOal ). are based 
on an analysis of the spatial correlation between ellipticities 
of galaxies. The study of weak lensing through the changes in 
magnitudes and apparent position in the sk y that is causes 
is an even more challeng ing measurement (|Scranton et al.l 
l2005l : IZhang fc Pen|[2005h . 

Most of the difficulties in measuring the weak lensing 
signal come from observational systematics such as uncer- 
tainties in the determinati on of the po i nt-spread function 
(PSF) and selec tion biases l|Kaiseij|200(il : lErben et al.ll200ll : 
I Vale et al.ll2004 ). There also exist astrophysical errors re- 
l ated to uncertainties in pho tometric redshift calibration 
(ilshak fc Hi rata 20051: IVan Wa erbeke et al.. ,2006, ') , intrinsic 
alignments due to physi cally associated lens-source pairs 
llMandelbau m et al.l l2006l ') or even distant lens-source pairs 
iHirata fc S cliak (20(3) • 

Most of the methods used to assess 
the importance of all these effects are theoretical or numeri- 
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cal. For the latter, one should also account, in principle, for 
uncertainties associated with numerical errors inherent to 
N-body simulations, neglected baryonic cooling effects and 
to a larger extent a poor understanding of the theory of 
galaxy formation. 

The most realistic simulations of weak gravitational 
lensing performed so far rely on a dark matter distribu- 
tion taken from an N-body simulation and galaxies drawn 
from a Halo Oc cupation Model (HOD) parameterized by 
halo mass alone (|Hevmans et alJbood : IVan Waerbeke et all 
I2OO6I ). without any information about apparent magnitudes 
or galaxy colors. Even if this information was available from 
such models, the mere fact that properties of dark matter 
halos depend not only on their masses, but also on their 
assembly histories — and t hat these latter have signific ant 
effects on galaxy clustering (|Croton. Gao fc WhitduOOfih — 
rules out HOD models based on halo mass as precision tools 
to model effects that rely heavily on a realistic description 
of the correlations between galaxy and dark matter distri- 
butions. The consensus is that an analysis which includes 
more realistic galaxy populations at various wavelengths are 
needed. 

This paper presents an approach which is a first step to 
fill this gap and perform more realistic weak lensing simu- 
lations thanks to a more realistic description of galaxy pop- 
ulations. It consists in building a tool which we call Lensed 
Mock Map Facility (lemomaf) hereafter, to link together 
three well established numerical techniques, each one tack- 
ling a different aspect of the problem. More specifically, a 
state-of-the-art semi-analytic model (SAM) tracks the prop- 
erties of galaxies within a high resolution N-body simulation 
as they evolve in time. Then a light cone is assembled from 
the outputs of the simulations to convert "theoretical" quan- 
tities into observables. Finally, a ray tracing algorithm ex- 
tracts the weak lensing signal from the cone. The SAM of 
galaxy formation that we use in th is work is GALICS (G ALax- 
ies In Cosmological Simulations, (jHatton et al.|[2003l )). and 
mock galaxy maps along with dark matter cones are ob- 
tained through a random tiling technique of simulat ion snap- 
shots , using the momaf (Mock Map Facility, (Blai zot et al.l 
I2OO5D ) pipeline. 

This paper is organised as follows. In Sec. [5] we review 
the characteristics of GALICS and momaf which are relevant 
to the present study. In Sec. [3] we explain the weak lensing 
formalism and its implementation in lemomaf. In Sec. |3] 
we present the statistics of the weak lensing convergence 
and the effect of weak lensing on differential galaxy counts. 
Finally, we discuss our results and outline prospectives in 
Sec. [5] 
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are turned into mock catalogues using MOMAF. The following 
section summarizes how this is achieved. 



2.2 MOMAF 



MOMAF (jBlaizot et al.|[2005l ) is a tool which converts theoret- 
ical outputs of hierarchical galaxy formation models into cat- 
alogues of virtual observations. The general principle is sim- 
ple: mock light cones are generated using semi-analytically 
post-processed snapshots of cosmological N-body simula- 
tions. These cones can then be projected to synthesise mock 
sky images. 

MOMAF u ses a random tiling technique described in 
iBlaizot et all |200i) to build mock observations from the 
redshift outputs of G alios. As explained in iBlaizot et al.l 
(|2005l ). several different observing cones can be generated 
from the same set of outputs of GALICS, by changing ei- 
ther the line of sig ht or the seed for the random tiling. 
IBlaizot et al.l (|2005l ) also discuss the bias induced by this 
random tiling approach on the clustering signal of the final 
maps. 

In this paper, we build 25 cones of galaxies and 25 cor- 
responding cones of dark matter with seeds and lines-of- 
sight randomly chosen for each observer position, using the 
same seed for every matching pair of galaxy and dark mat- 
ter cone. These 25 cones allow us to infer an estimate of 
the dispersion of clustering measurements, that is to say, an 
estimate of the cosmic variance associated with our mock 
catalogues. We use galaxy magnitudes in two filters: SDSS 
r and JOHNSON K, because the results in these ban ds have 
alre ady been thor o ughly examined by IBlaizot et al. I (l2005l ) 
and IBlaizot et al.l (I2OO6I ) . The results of the galaxy counts 
from GALICS-MOMAF and its comparison to observations are 
shown in Fig. [T] and galaxy redshift distributions for dif- 
ferent cuts in magnitude are shown in Fig. (2] In order to 
mak e semi-analytic predictions of the weak lensing statis- 
tics l|Van Waerbeke et al.|[200lh for this broad distribution 
we describe it with the following functional form: 



exp 



(1) 



where r(a;) is the Gamma function and a, (5 and zq are free 
parameters to determine for each distribution. 

However, we note that given the rather small size of the 
original N-body simulation box (100 Mpc on a side) 
that was used to run GALICS on, our estimate of the cosmic 
variance is likely to be biased and has to be taken as a 
lower boundary on the true cosmic errors. Finally, the mock 
catalogues we will use in the following are 4400 h-^ Mpc m 
depth, and 1.4° x 1.4° in angular size. 



2.1 GALICS 

GALICS is a hybrid model of galaxy formation which com- 
bines cosmological dark matter N-body simulations with a 
semi-analytic descr iption of bary onic processes. The model 
is fully described in lHatton et al.i (2003), and the version we 
use here is the sa me as that used in the previous papers of 
the GALICS series ( Hatton et "al]|2003l : iDevriendt et al.|[2005l : 
IBlaizot et al. I I2OO4I ). We briefiv recall what are the main in- 
gredients in Appendix 13 and |B] Eventually, GALICS outputs 



3 LENSED MOMAF 

3.1 Weak Lensing Formalism 

In this section we provide a summary of the weak lens- 
ing equations relevant for building lemomaf an d refer the 
reader interested in a more detailed treatment to I Jain et al.l 
(l2000l) . 

We use comoving coordinates and place ourselves in the 
framework of an isotropic and homogeneous universe that 
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Figure 1. Galaxy co unts in tlie filters SDSS r l lYasuda et alj I2OOII ) and JOHNSON K l lDiorgovski et alj Il995l : [Gardner et alj Il996l : 
iMoustakas et aLlll997^ obtained through the GALICS-MOMAF pipeline. Squares are the observational points, and lines the results of the 
analysis of several virtual light cones built with the simulation. On the right hand side panel one can clearly see when the incompleteness 
in the number of low luminosity mock galaxies kicks in due to the finite mass resolution of the simulations. 




Figure 2. Distributions in redshift for different cuts in r and K magnitude of the galaxy cones. The dots show the mean value of ns{z) 
obtained by averaging over the 25 cones, and the error bars show the 1-ct error on this mean value. The continuous line is the fit to the 
function in Eq. [T] 



can be described by the Robertson- Walker metric, where in 
the presence of a perturbative gravitational potential the 
change in a photon's direction can be written as: 



located at coordinate perpendicular to the line of sight, 
a deflection of 



da = -2V±<^dx (2) 

where da is the photon's deviation, is the gravita- 
tional potential, is the gradient in the direction perpen- 
dicular to the photon line of propagation and x is the radial 
comoving coordinate. A deflection at x' produces in a plane 



dx = r{x-x')daix') (3) 

where r(x) is the comoving angular distance, that in a plane 
universe (^a + = 1) is equal to x- Integrating along the 
perturbed trajectory of the photon, and then dividing by 
r(x) we obtain the angular position at position x 
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Six) = ^"(0) r dx'rix - x')Vx<^ (4) 

nx) Jo 

This treatment holds for a single photon. To obtain the 
effect on an extended object we calculate the Jacobian of the 
transformation {ddi{x)/d6j{0)), which we call Aij where i 
and j denote perpendicular directions in the plane perpen- 
dicular to the line of sight, in which 9 is also located: 

Jo '■(X) 
the Jacobian matrix is usually decomposed as follows: 

V —72 + 1 — K + 7iy 

where n is the convergence, 7 is the shear and to the rotation. 

The relation between the gravitational potential and 
the mass density perturbation 5 = p/p — 1 is given by: 

where Hq is the Hubble constant at the present, c is the 
speed of light and a is the expansion factor. 



3.2 Multiple plane formalism 

The multiple plane formalism consists in dividing the space 
between the source and the observers into A'^ equidistant (in 
comoving coordinates) planes perpendicular to the l ine o f 
sight. Following the convention of liamana & Mcllicj (1200 ll ) 
we place ourselves in a cartesian coordinate system noted by 
{xi,X2,y) where the origin is the observer and y indicates 
the direction of observation. For small angular size fields of 
view, y can be identified with the radial comoving distance 
X- The inter-plane comoving distance will be called Ay. The 
projected density contrast on the i-th plane located at yi is 
given by: 

rvi 

T,i{xi,X2)= / dyS(xi,X2,y) (8) 

Vi-l 

which defines a corresponding two dimensional potential 
^' = 2 J dycj). The position of a light ray in the n-th plane 
can then be found using the multiple plan lens equation: 

e„ = k - V (9) 
^^^j_gKx:,^xOiMu,A. (10) 

where I is the identity matrix and Ui is defined by: 

/an*. 9i2*A , . 

" \^2l■i'^ 922*J ^^^^ 

where dij symbols stand for partial differentiation. 



3.3 Ray Tracing Algorithm 

The algorithm can be divided in three steps. First, the pro- 
jection of the dark matter distribution onto the planes. Sec- 
ond, the calculation of the potential and its first and second 
derivatives in these planes and finally the ray tracing from 
the observer to the last plane. In practice, things proceed as 
follows: 

1. Once we have the dark matter particles of the sim- 
ulation positioned inside an light cone build with momaf, 
we define an orthogonal coordinate system {xi,X2,y) at the 
origin of the cone, where the y axis is directed from the ob- 
server along the symmetry axis of the cone. We then project 
(along the y axis) all the particles onto A'^ equidistant planes 
parallel to the one located at the origin of the cone, which 
means for example that particles originally between (xi, X2, 
y) and (x\,X2,y 4- Aj/) now sit at {x\,X2,y -\- Ay). 

2. In each plane we interpolate the overdensity [5 in 
equation 0), using a cloud-in-cell (CIC) algorithm, over an 
uniform square grid, G, of Ng cells on a side of comoving size 
Lside- Ng and Lside are kept identical for all planes. In order 
to solve the Poisson equation we pad with zeroes around the 
grid G, thus creating a grid of size 2 x Ng in side. We then 
perform a fast Fourier transform (FFT) of the overdensity 
grid, and in Fourier space we divide by —l/k^ to inverse 
transform and obtain the gravitational potential 'i! on the 
grid. Finally, using finite differences methods we obtain the 
first and second derivatives of the potential along the xi 
and X2 directions. This allows us to construct the matrix U 
in Eq. [TT] for each point in the plane. A det ailed numerical 
implem ent at ion for this step can be found in lPremadi et al] 
(|l998h . 

3. We initialize the ray positions over a uniform grid on 
the nearest plane to the observer. For each ray we interpolate 
the values of a and the elements of U using their tabulated 
values on the CIC grid. Applying equation ([9)l we then figure 
out the position of each ray on the next plane, and again 
interpolate the values of a and U at this new position to 
calculate A, and so on and so forth until we reach the final 
plane. We store all these values, for each ray on every plane. 

3.4 Shearing of galaxies 

The galaxies in GALICS are represented by three components: 
disk, bulge and burst. Geometrically speaking, on a mock 
image, the disk is seen as an ellipse, and the bulge and the 
burst as circles. The burst is treated as a punctual source 
and only its magnitude will be modified by weak lensing. 
The circle and the ellipse can both be parametrised by their 
shape matrix, according to their weighted quadrupole mo- 
ments. 

For an ellipse of major axis a, minor axis h and with the 
major axis making an angle /3 with respect to an horizontal 
reference axis, the shape matrix can be written as: 

_1_ /a^ cos^ /? + 6^ sin^ /3 (a^ - 6^) sin/3cos /3 \ 
^ ^ 7ra6 V(a2 - 62)sin/3cos/3 sin^ /3 -I- fe^ cos^ /3 ) 

(12) 

The case of the circle is recovered by setting a — b and 
/3 = 0. Supposing that the Jacobian matrix does not vary 
much along the corresponding surface of the galaxy (which 
is a good approximation in the case of weak lensing caused 
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by large scale structure) the lensed shape matrix Q' can be 
written as: 



Q' 



(13) 



In this way the lensed properties of the disk and the 
bulge can be easily found if one knows the Jacobian matrix. 
In LEMOMAF the A matrix used to lens a galaxy is the av- 
erage of four matrices, located at the point of impact of the 
four closest rays on the nearest plane to this galaxy. We lens 
bulge and disks separately. 

This is used to produce mock lensed images useful for 
a more realistic treatment of the detection in simulations of 
the lensing signal through cosmic shear. In this paper we 
won't use this capability. We will infer the measurements 
of K directly from the numerical values used to modify the 
galaxy's properties. 



3.5 Limits of the method 

Most of the limitations of our ray tracing implementation 
are resolution issues that come from the use of N-body sim- 
ulations and the interpolation grids used to trace the light 
rays. These errors in the calculation of the weak lensing sig- 
nal from numerical simulations have been thoroughly quan- 
tified before ( Jain et al, 2000 ; Valefc White 2003). We recaU 
these results here, explicitly pointing out the place where 
the parameters that we use limit the approach the most. 

Two relevant scales can be defined to asses the quality of 
the weak lensing signal in the simulation. The first one, Ug, 
relates to the size of the Fourier grid where the lensing signal 
is obtained, and the second, (t„, to the finite resolution of 
the N-body simulation. We define these quantities as: Gg — 



jNg 



JK 



1/3 



where L, 



is the size of the 



grid in the plane where we interpolate, 
cells in one dimension, Li,ox is the size of the simulation box 
and Npart is the number of particles inside this box. 

If (jg is larger than (t„, the power of the signal on small 
scales will be dominated by the Fourier grid cut-off. If it is 
smaller, the power of the signal is dominated by noise in the 
N-body simulation. Moreover, if ag is not only larger than 
CT„ but also larger than any significant physical scale in the 
simulation, the features of the overdensity field that produce 
the deflection of the rays on that scale are wiped out, and no 
weak lensing signal is measured by the ray tracing method. 

We choose a value of ag/an ~ 1, guided by the results of 
previous studies. In our case, since we use a simulation with 
Lbox = 100 Mpc and A^p — 256'', this translates into a 
resolution of (7„ — 390 kpc. Our ray tracing simulation 
uses a grid with Lside ~ 200 Mpc and Ng — 800, which 
translates into ag = 250 h~^kpc. 

For a source at a given redshift, the weak lensing signal 
probes structures at intermediate redshifts between source 
and observer. This can be seen from the equation (|10p where 
the distance combination r(xn — xO^(xO/''(Xn) plays the 
role of an efficiency function for the lensing convergence. 
This distance combination peaks at an intermediate redshift 
between the observer and the source as shown in the upper 
panel of Fig. [3] where it is plotted normalized to the plane 
distance for three different source positions. 

When one takes into account the redshift distribution of 
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Figure 3. Top Panel: Distance combination in equation [To] nor- 
malized to the distance to the source. The three curves correspond 
to a distance to the source of 1000, 2000 and 3000 /t~^ Mpc. Mid- 
dle panel: Sum over the curves of top panel weighted by ns{z) at 
the redshift of the source plane. The solid curve, labeled as BASE, 
shows this sum for a population of galaxies homogeneously dis- 
tributed in redshift. The curves labeled as JDHNSDN K and SDSS 
r show the sum for the respective galaxy distribution in Fig[2] 
Lower panel: value of the angular resolution in our ray tracing 
simulation as a function of distance from the observer. 



the sources used to measure the weak lensing signal, ns{z), 
one realizes that each redshift contributes a different amount 
to each plane. In order to estimate the joint contribution 
from the lens efficiency and the redshift distribution of the 
sources, we consider that each efficiency curve is multiplied 
by the value of ns{zi) where Zi is the redshift at which the 
efficiency curve is calculated, and then we add all the differ- 
ent efficiency curves together. In the middle panel of Fig. [3] 
we show the results of this calculation for a uniform ns{z) 
as well as for the six source distributions presented in Fig. 


From the middle panel of Fig.|3] one can see that the dis- 
tances probed by galaxy populations measured in the SDSS 
r filter are comprised between — 1000 Mpc (redshifts 
to 0.5), while for the JOHNSON K filter this distance interval 
lies between 1000-3000 h'^Mpc (redshifts from 0.5 to 2.0), 
with a broad peak around 1500 h~^Mpc. 

Once the value of Og is fixed, we can determine the 
angular resolution as a function of the distance from the 



6 J.E. Forero- Romero et al. 



observer, as shown in the lower panel of Fig. |3]for our cho- 
sen value of (Tj. From this figure, one can check that in the 
redshift range probed by galaxies seen in the SDSS r filter 
this angular resolution is on the order of ^ 1.5 arcmin, and 
for that probed by galaxies in the JOHNSON K filter of about 
~ 0.5arcmin. 

Discontinuities induced by the random rotations and 
origin shifts of the tiled boxes in the construction of the 
cone, using planes that are not perpendicular to the path of 
the light rays prod uce artifacts that were also analyzed by 
IValefc White! (|2003l ). concluding that their effects on the lens- 
ing signal are negligible. Finally, other errors arise from the 
translation of dark halo mass resolution into completeness 
limits at a given magnitude in a given waveband at a given 
redshift for the galaxy population, or from reduced spatial 
correlations caused by the random tiling technique. All this 
limitations have been studied in d etail liatton ct al. ( 2 0031) . 
iBlaizot et all l|2004 ). [Blaizot et all l|2005f ) and lBlaizot et all 
1 20061 ) . In this paper we remain inside these limits to 
draw our conclusions, but explicitly indicate when they are 
reached. 



4 RESULTS 

In this study, we make two kinds of numerical experiments 
with LEMOMAF: (i) we measure the weak lensing signal in- 
duced by the dark matter cone at galaxy positions, and (ii) 
we study the change in galaxy counts caused by weak lens- 
ing effects. Before exploring the results of these experiments 
we briefly discuss the behaviour of the lemomaf ray tracing 
module. 



4.1 Ray Tracing Results 

In this section, we focus on the measurements of the mo- 
ments of convergence (k) to present the results of the ray- 
tracing module of lemomaf. 

The dark matter cones we use to make these measure- 
ments are 1.4° on a side, and have a depth of 4400 Mpc, 
with a distance of 100 h^^ Mpc between lensing planes. In 
each of these planes, we use a grid with size L^ide = 200 
Mpc which is split in Ng = 800 cells. Weak lensing proper- 
ties and galaxy counts are measured only in a central field 
of 1.0° on a side. We characterize k through moments of its 
distribution function such as its variance (ng), its skewness 



Ml 



and kurtosis 



Si{e) = 



(14) 



(15) 



We calculate these moments after smoothing the conver- 
gence maps with a circular top-hat filter of angular scale 
9. 

We are especially interested in two results to validate 
our ray tracing code: 

1) The dependence of (k^)^''^ and S3 as function of the red- 
shift of the source, for which analytic expressions can be 
derived. 




0.1 



1.0 



Figure 4. Convergence statistics as a function of redshift . Top 
panel: variance (divided by 10~'^). Bottom panel: skewness (di- 
vided by 10^). The solid lines are the numerical results. Each line 
represents a different value obtained after smoothing over a dif- 
ferent angular scale 6 regularly spaced every arcminute from 1 
to 9 arcminutes. The dotted line shows the theoretical trend pre- 
dicted by Bernardeu et al. 1997. The numerical trend only weakly 
depends on the smoothing scale 9. The numerical and theoretical 
results compare fairly well at ^ > 1, but at z < 1 the numerical 
results predict a logarithmic slope about two times steeper than 
that of analytical theory. 



2) The dependence of {^1)^'^, S3 and ^4 on the smoothing 
scale 0, for a given redshift of the source, for which published 
values exist in the literature. 

In Fig|3]we show our measurements for {lY)^^^ and S3 
as a function of redshift, for nine different smoothing angles 
9 spaced every arcminute from 1' to 9'. From this figure we 
can see that our computational values follow fairly well the 
expected theoretical trend. The agreement is better when 
the source plane is located at redshifts larger than z = 1. 
For closer source planes, located at 2; < 1, the slope of our 
numerical relation is steeper by about a factor of two. 

In Fig. [5]we plot our results for (k^)^''^, S3 and 5*4 as a 
function of the smoothing angle 9, for sources located in a 
single plane at z = 1 and a field of view of 1° x 1°. The order 
of magnitude of S3 and S4 is consisten t with a comp i lation 
of results for these moments made by IValefc White! (l2003t ) 
for smoothing scales of 4'. These authors quote an average 
of 5*3 ~ 135 ± 10 and S4 ~ (3.5 ± 0.5) 10* for measures ob- 
tained using different methods, with error bars reflecting the 
dispersion amongst reported values. For the same smooth- 
ing scale we find S3 = 115 ± 5, and 54 = (2.5 ± 0.3) 10*, 
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Figure 5. Convergence statistics as a function of smoothing scale. 
Sources are all located in a single plane at z = 1. Top panel: 
variance (divided by 10~^). Middle panel: skewness (divided by 
10'^). Bottom panel: kurtosis (divided by 10"*). Error bars indicate 
the 1-iT dispersion measured around the mean for the 25 cone 
realisations. 



keeping in mind that uncertainties coming from the numer- 
ical simulations themselves can contribute for an extra 10% 
error in our c ase, a nd that theoretical methods compiled in 
IValefcWhitel (|2003l ) suffer from a similar underestimate. 

Being reasonably confident that we are able to produce 
reliable weak lensing measurements, we can now couple our 
ray tracing code to the outputs of a hybrid model of galaxy 
formation to study the source-lens clustering efi'ect on the 
convergence statistics, and the effect of weak lensing on 
galaxy counts. 



4.2 Source Lens Clustering 

In the previous section we measured the statistics of the 
convergence of rays that uniformly covered a fraction of the 
sky. In reality, we only have access to the convergence signal 
measured for galaxies that have specific clustering proper- 
ties and, even more important, are correlated to the lensing 



potential. This is known in the literature as the Source-Lens 

Clustering (SLC) effect. 

Theoretical studies of the SLC effect (|Bernardeaulll998h 
predict that it alters the higher order statistics of the con- 
vergence. More precisely, while (k^)^^^ should remain un- 
changed according to theory, S3 should have a lower value. 
As a matter of fact, S3 is known to be sensitive to al- 
most independently of erg, and a combined analysis of the 
skewness and the variance of the convergence could in prin- 
ciple provide new constraints on the values of Qm and as- 
Hence the cosmological interest to accurately quantify the 
inffuence of the SLC effect on the higher order statistics of 

K. 

Previous numerical work on the SLC effect 
(iHamana et al. has focused on its impact on S3 

estimations. This work was carried out using a simple bias 
model to paint a population of galaxies on top of the dark 
matter density field. To our knowledge, our work is the 
first attempt to use a galaxy population self-consistently 
derived from a N-body simulation to estimate of the impact 
of the SLC effect. More specifically, with momaf, we obtain 
both a dark matter distribution and its matching galaxy 
distribution in a light cone, this latter being derived by 
post-processing the dark matter simulation with the GALICS 
semi-analytic model. We then calculate the weak lensing 
signal over all the field at all the planes used in the ray 
tracing simulation with lemomaf, and use this information 
to shear the properties of galaxies which are in the cone, 
while storing the value of the convergence computed at 
each galaxy position. In this section, we therefore return 
to the analysis of the convergence statistics, but only for 
the subset of values measured at each galaxy position by 
interpolation of the values at the neighboring rays. 

To quantify the impact of using a self-consistent mod- 
elling of the galaxy population on our results, we couple 
the dark matter and galaxy cones in three different ways. 
In the first way, that we call MATCH, we shear the galaxies 
according to their true underlying dark matter distribution. 
In the second way, called RANDOM, we keep the same pair 
of galaxy cone/DM cone as in the MATCH case, but erase 
some of the spatial clustering information by randomizing 
the positions of galaxies over the sky plane while keeping 
their distance to the observer constant (i.e. we preserve the 
galaxy redshift distribution). Finally, in the third case, called 
NO MATCH, we use the full spatial clustering information for 
galaxies, but shear them using a different underlying dark 
matter distribution from the one with which their proper- 
ties were derived. The idea behind these "tricks" is that 
cross-comparisons between the MATCH, NO MATCH and RANDOM 
methods should provide us with a better understanding of 
where the impact of the source-lens clustering effect on the 
convergence statistics comes from. 

We use two broad band luminosities for each galaxy, 
those measured in the SDSS r , and JOHNSON K filters. We 
build 25 different light cones in order to minimize the bias 
effects induced by the random tili ng technique, and m axi- 
mize the accuracy of the statistics (iBlaizot et al.ll2005l ). For 
each light cone we output both the gala:xies and the dark 
matter distribution. To mimic observational effects as best 
as we can, we select galaxies on which the lensing signal is 
to be measured based on their apparent magnitude. The re- 
sulting redshift distributions are shown in Fig. [2] where we 
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Figure 6. Comparison of the convergence statistics for three different magnitude ranges in the SDSS r filter, and for three different 
methods of measuring the convergence at each galaxy position. The solid line indicates results obtained with the RAND method, the dashed 
line with the NO MATCH method, and the dotted line with the MATCH method, {k?) has been divided by 10"'*, ^4 by 10*. Overplotted 
diamonds show values computed with the semi-analytic model described in Van Waerbeke et al. 2001. 
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Figure 7. The R{S) factor as defined in Eg. llOl for {k^) and S3, for the three different magnitude ranges in the SDSS r filter. The solid 
line compares the methods MATCH and RAND. The dashed line compares the methods MATCH and NO MATCH. 
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have arbitrarily split magnitudes in three bins: 20 < m < 21, 
21 < m < 22 and 20 < m < 22. Once again, we recall that 
each light cone is 1.4° on a side, and has a depth of 4400 

Mpc, with a distance between lensing planes of 100 
Mpc. In each of these planes we use a grid of size Lside = 200 

Mpc split into Ng = 800 cells. Weak lensing properties 
and galaxy counts are measured in a centered field of 1.0° 
on a side. 

Fig. |6] and |8] show the results for the higher order mo- 
ments of the convergence and the six redshift distributions 
plotted in Fig. [2] Each linetype corresponds to a differ- 
ent method to construct the maps; RANDOM, NO MATCH and 
MATCH. We also compar e our numerical results to th e semi- 
analytic calculations of IVan Waerbeke et al.l (j200ll 'l. Over- 
all, the semi-analytical predictions of these authors show 
very good agreement with the lensing signal measured in 
our mock catalogues. Moreover, the better agreement with 
the RANDOM mocks is somewhat expected: out of the three 
methods we presented here, this is the one which follows the 
most closely the assumptions made in IVan Waerbeke et al.l 
l|200ll l. Most of the small scale deviations from the semi- 
analytic trend in the variance plots can be attributed to the 
finite spatial resolution of our simulation (see Section 3.5). 

We introduce the parameter R as done by 
iHamana et al.l (|2002l ) to quantify the amplitude of the 
SLC effects: 



4.3 Galaxy counts 

Weak lensing enlarges the area of the sky which is ob- 
served, thus lowering the number density of objects which 
are detected. In the same time, it causes galaxies to ap- 
pear brighter, thus increasing their number density for a 
given apparent magnitude. The net effect depends on the 
slope s(m) of the number counts of galaxies. Let us write 
No{m)dm as the number of galaxies with magnitudes be- 
tween m and m + dm, and s(m) = d\og No(m) / dm. If the 
sources undergo a magnification ji = 1/ det A (see section 
3.1 for a definition of the matrix A), and N'{m) is the mag- 
nified number of sources corr esponding to No(m), we can 
write Af'(m) = ^.^-^"-^ Npi m) l|Broadhusrt et. al.lll995l : |jainl 
l2002l : IScranton et aklbOOs h . 

The effect of the magnification on the solid angle, which 
is responsible for a factor can only be detected if we use 
a grid that allows us to have a resolution in the order of 0.1 
arcmin to make a proper estimation of the deflection angle 
for a given galaxy. The angular resolution of our simulations 
is roughly 1 arcmin, which is clearly not enough. Therefore 
we expect that only the lensing effects on galaxy magnitudes 
(which does not require high angular resolution to be seen) 
will be perceptible. This simplifies the expression of mag- 
nified number of sources to N'{m) = fj,^'^" No{m), which in 
turn reduces, in the weak lensing regime {fi = 1 + 2k), to: 



R{S) = 



5, 



(16) 



N'{m) 
Noim) 



1 = 5.0 K s 



(17) 



where S can be (k^) or S3, and NO SLC is one of the methods 
RAND or NO MATCH. Fig. [7] and |9] show the results for this 
expression. 

From these figures, it is clear that the results for {k^}, 
Ss and S4 in both filters, and in the three different mag- 
nitude bins show the same trend: these statistics are sys- 
tematically lower for the MATCH case than the for the NO 
MATCH and RANDOM cases, with these latter being remarkably 
similar. Furthermore, we clearly see that the effect is more 
pronounced for the narrower redshift distribution of sources 
(i.e. SDSS r band). These results, which are obtained using 
a fairly realistic galaxy distribution, are in broad agreement 
wit h the simpler approach advocated bv lBernardeaul (|l998l ) 
and iHamana et all (j2002h . However, we also find a signifi- 
cant SLC effect on the two point statistics {R ~ 20% for 
the SDSS r distributions, R ~ 5% for the JOHNSON K dis- 
tributions), which is not expected from perturbation the- 
ory alone. We suspect that it arises because our simulation 
probes the highly non-linear regime. A comparison of the 
MATCH to the NO MATCH run — where the clustering of sources 
is preserved albeit galaxies positions are not correlated with 
that of the underlying dark matter distribtion — strongly 
suggests that the signal comes from the adequation of the 
intra halo galaxy population with the depth of its host po- 
tential well. 

Finally, we would like to emphasize that such a mock 
catalogue approach should allow a fast calculation of the 
SLC effect for future weak lensing surveys intending to use 
cosmic shear as a precision cosmology tool. 



The numerical experiments we carry out with lemomaf 
are meant to explore the effect of weak lensing on galaxy 
counts via the estimate of this ratio, as we vary the size of the 
field in which the measurement of the counts are performed. 

In the 25 original fields of 1° x 1° we compute galaxy 
counts both for unlensed and lensed fields (using the MATCH 
and NOMATCH methods) . We then cut each of these fields into 
smaller square patches of angular size 12', 15', 20', 30' and 
1° on a side, and measure the counts in all of these sub-fields. 
We also calculate for each patch the ratio N{m) /No{m), and 
organize our results as follows. For a patch of size 6, taken 
from a field of original size O we have A'p = (O/Q)^ patches. 
Labeling subfields as ^ = 1, . . . , A'p in every original uncut 
O field, we then proceed to stack together the 25 patches 
which have the same index £. For each patch in the stack, 
we estimate the ratio of lensed to unlensed integrated counts, 
and compute the mean for each stack. 

In Fig. 1101 and Fig. [Tl]we plot these means for our Np 
stacks, for each value of 9 (column panels) and for both 
methods, NO MATCH and MATCH (top and middle row panels 
respectively). In the bottom row of these figures, we show 
A = 5.0 s K, in order to facilitate the interpretation of the 
behaviour of the lensed to unlensed counts ratio in terms 
of Eq. 1171 In this intent, we estimate k using an average 
value measured for the galaxies which lie in the magnitude 
bin of interest in the uncut field of size O. From our ray 
tracing results (see Fig. |4]) we know that we are in the weak 
lensing regime where k ^ 1, so Eq. [T7] tells us that we can 
expect a depletion in galaxy counts only when (k) and s 
have different signs. 

In the SDSS r band (Fig. llOp . A (in percent) is restricted 
to the interval [0.0, 1.0]. High fluctuations in the counts ratio 
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Figure 10. Percentile difference between lensed (A'^') and unlensed (A^o) galaxy counts in the SDSS r filter, calculated with the MATCH 
(upper panel) and NO MATCH (middle panel) methods. The lower panel shows the values of A = 5.0 s {Kgai)hm^ where s is the 
logarithmic slope of the counts and {Kgai)bin is the average value of k measured from the galaxies in that magnitude bin over the field 
of side size l°.Each vertical line splits the plot in panels which contain results for different patches of angular size 9 on a side, cut inside 
an original field of size G = 60' on a side. Thus the number of curves in each row of panels is Np = (0/9)^, and each curve represents 
the mean of the measurement over 25 uncorrelatcd patches. The rightmost upper and middle panel also show the l-a dispersion between 
the 25 © fields [Np = 1) as error bars. 




Figure 11. Same as Fig. llOl for the JOHNSON K filter. The behaviour for magnitudes A" > 21 is a numerical resolution artifact which leads 
to an incompleteness in the number of low luminosity galaxies (see also right panel in Fig[TJ 



for the low values of the magnitude are due to an ever smaller 
number of bright galaxies being present in a subfield when 
the subfield size decreases or the brightness of the source 
increases. In the magnitude interval rrir — [19, 22] where 
this efi'ect becomes negligible, we see that galaxy number 
counts can be enhanced by up to 1%. 

In the JOHNSON K fiher (Fig. [TTJ, A (in percent) takes 
values in the interval [—4.0, 4.0], although its change of sign 
between the magnitudes 21 and 23 is entirely due to the 
mass resolution of our N-body simulation which translates 
into an incompleteness for galaxies fainter than tuk = 21 
(see also Fig. [T} . We therefore restrict ourselves to the mag- 
nitude interval rriK ~ [16, 21] for which there exists a good 
agreement between modeled and observed counts in terms 



of slope, and the number of bright galaxies per subfield is 
high enough. In this magnitude interval, the counts can be 
enhanced up to a 3%, for high values of niK- 

In both filters we note that faint galaxies are always 
more enhanced than bright ones: A as a function of a mag- 
nitude is a monotonically increasing function. Of course, this 
trend has to break down at some point, when the slope of the 
faint galaxy counts turns over, as shown in Fig. 1111 even if in 
this case it is purely an artifact due to finite mass resolution 
in our N-body simulation. 

Unsurprisingly, the dispersion around the theoretically 
expected ratio increases when the angular size of the field 
decreases. Finally, Fig. llOl and llll show that the source lens 
clustering effect does not play an significant role in enhanc- 
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ing galaxy number counts: the NO MATCH and MATCH methods 
pretty much yield the same quantitative results. 

However, from this experiment we confirm our suspicion 
that the modification of the solid angle is not resolved in our 
simulations, consequently higher angular resolution must be 
attained if one hopes to use lemomaf in the study of angular 
correlations induced by cosmic magnification. 



5 DISCUSSION AND CONCLUDING 
REMARKS 

In this paper we presented the Lensed Mock Map Facility 
that combines dark matter N-body simulations and an hy- 
brid model of hierarchical galaxy formation to make mock 
lensed images and convergence maps, thanks to a ray tracing 
algorithm. More specifically, the results presented here were 
obtained using a cosmological N-body simulation performed 
with a treecode, the GALICS model of galaxy formation, the 
MOMAF pipeline which constructs galaxy and dark matter 
light cones, and a ray tracing algorithm through multiple 
planes to account for the weak lensing effect. This tool suf- 
fers from all the shortcomings inherent to each of these tech- 
niques. However, for GALICS and momaf these limitations 
have been carefully identified in a series of papers published 
over the past four years. As far as the ray tracing algorithm 
is concerned, its limitations in terms of the N-body simula- 
tion parameters have also been discussed quite thoroughly 
in the literature. To sum up, these shortcomings impose a 
limitation on the size of the fields that can be constructed, 
and the angular resolution that can be reached in the con- 
struction of the lensed images and the convergence maps. 

Working within the proper interval of validity of these 
methods, we performed two numerical experiments with 
LEMOMAF. The first one measured the convergence signal 
induced by the dark matter density field at galaxy posi- 
tions in a light cone. Different methodologies for this mea- 
surement were implemented, in the aim of testing the con- 
sequences of the source-clustering effect on the probabil- 
ity density function of the convergence. We found that the 
SLC effect skews the PDF towards lower values of the con- 
vergence and that, in some cases, it makes this PDF look 
more gaussian than that obtained without including SLC, 
as expected from theoretical considerations. However, even 
when probing the same dark matter distribution, the pre- 
cise trend of the SLC effect depends sensitively on the spe- 
cific distribution of the galaxies we consider. For instance, 
we demonstrated that a narrower redshift distribution of 
the sources is more sensitive to the SLC effect. This could 
be problematic for future lensing surveys which intend to 
perform shear measure ment in thi n redshift slices, a tech- 
nique called tomoaravhv IHi] (|l999l ). For the JOHNSON K 
filter, the SLC effect has an impact at the few percents level 
(2 — 5%) on the estimations of ag from two point statistics. 
Thi s level of contamination was neglected i n previous analy- 
sis l|Bernardeaull 19981 : Ih amana et al.ll2002l ). because its am- 
plitude was well below the lar gest we ak lensing surveys ac- 
curac y at that time (yiRMO S: IVan Wa orbcko ct al. (20q3), 
RCS: iHoekstra et~aLl |20o3)). However, future -nearly full 
sky- missions like SNAP and LSST will have to reach a pre- 
cision of 10""^ on the shear measurem ent, i.e. roughly 0.1% 
from on the shear two-points statistics l|Van Waerbeke et al.l 



(^ 20061 )). At this level of precision, the SLC will be a major 
source of systematics, and the only way to tackle this issue 
is to have photometric redshifts for each individual galaxy, 
sources and lenses. This strengthen the requirement that 
future lensing surveys will have to cover a wide range of 
the optical spectrum, from U band to near infrared, with 
narrow band filters, similar to the COMBO-17 approach 
iHevmans et al.ll2004l ). 

The second numerical experiment measured the lensed 
to unlensed galaxy counts ratio. The value of this ratio was 
obtained for various angular sizes of observational fields. The 
general trend of the results in the simulation can be under- 
stood using Eg. 1171 where only magnitude changes played an 
important role in enhancing the counts in our simulations. 
We learnt from this experiment that in order to have a real- 
istic treatment of the magnification effects over the change 
of galaxy positions in the sky we need to go up in angular 
resolution. 

Among the ideas that remain to be investigated using 
LEMOMAF are those that take advantage of mock images at 
multiple wavelengths to identify the best strategies for mea- 
suring the shear, as well as those which intend to study the 
bias introduced by intrinsic alignments or other systematic 
effects on this measurement. Future prospects with LEMO- 
MAF include the simulation of galaxy-galaxy lensing and cos- 
mic magnification. This kind of signal demands a resolution 
in our simulations in the order of 0.1 arcmin. With the N- 
body simulation we used in this paper, and FFT methods 
we barely achieve a resolution of ~ 1 — 2 arcmin. In order to 
reach higher resolutions we do not plan to only rely on more 
resolved N-body simulations but also to switch ray tracing 
strategy and look in the direction of smooth particle lensing 
(|Aubert et aLlbOOd ). as we have gained confi dence from re- 
cent studies with GALICS (|Blaizot et aklbOO^ ) that the small 
scale distribution of galaxies produced from an adequate res- 
olution N-body simulation and a new positioning scheme of 
galaxies inside the halos can be accurate enough to attain 
this goal. We also plan to use lemomaf to help design future 
lensing surveys, which will need to employ a tool taylored 
to tackle non-linear effects such as SLC. 
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APPENDIX A: DARK MATTER 

The cosmological N-body simulation (|Nininl [l999l ) we use 
throughout this paper assumes a flat Cold Dark Mat- 
ter cosmology with a cosmological constant {Qm ~ 1/3, 
JIa ~ 2/3), and a Hubble parameter h = _ffo/[100 km s~^ 
Mpc~^] = 0.667. The initial power spectrum was taken 
to be a scal e-free (rta = 1 ) pow er spectrum evolved as 
predicted by I Bardeen et al.l (|l986l ) and normalised to the 
prcsont-day abundance of rich clusters with erg = 0.88 
( Eke, C ole, & Frenk 1996 ) . The simulated volume is a cubic 
box of side Lb = 100/i~^Mpc, which contains 256'^ parti- 
cles, resulting in a particle mass rrip = 8.272 X 10^ Mq and a 
smoothing length of 29.29 kpc. The density field was evolved 
from z — 35.59 to present day, and we out-putted about 100 
snapshots spaced logarithmically with the expansion factor. 

In each snapshot, w e identify halos u sing a friend-of- 
friend (FOE) algorithm l|Davis et al.lll985^ with a linking 
length parameter b = 0.2, only keeping groups with more 
than 20 particles. At this point, we define the mass Mfof 
of the group as the sum of the masses of the linked parti- 
cles, and the radius Rfof as the maximum distance of a 
constituent particle to the centre of mass of the group. We 
then fit a tri-axial ellipsoid to each halo, and check that the 
virial theorem is satisfied within this ellipsoid. If not, we 
decrement its volume until we reach an inner virialised re- 
gion. From the volume of this largest ellipsoidal virialised re- 
gion, we define the virial radius R^ir and mass M^ir . These 
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virial quantities are the ones we use to compute the cool- 
ing of the hot baryonic component. Once all the halos are 
identified and characterised, we build their merger history 
trees following all the constituent particles from snapshot to 
snapshot. 



APPENDIX B: LIGHTING UP HALOS 

The fate of baryons within the halo merger trees found above 
is decided according to a series of prescriptions which are 
either theoretically or phenomenologically motivated. The 
guideline - which is similar to other SAMs - is the following. 
Gas is shock-heated to the virial temperature when captured 
in a halo's potential well. It can then radiatively cool onto a 
rotationally supported disc, at the centre of the halo. Cold 
gas is turned into stars at a rate which depends on the dy- 
namical properties of the disc. Stars then evolve, releasing 
both metals and energy into the interstellar medium (ISM), 
and in some cases blowing part of the ISM away back into 
the halo's hot phase. When haloes merge, the galaxies they 
harbour are gathered into the same potential well, and they 
may in turn merge together, either due to fortuitous col- 
lisions or to dynamical friction. When two galaxies merge, 
a "new" galaxy is formed, the morphological and dynami- 
cal properties of which depend on those of its progenitors. 
Typically, a merger between equal mass galaxies will give 
birth to an ellipsoidal galaxy, whereas a merger of a massive 
galaxy with a small galaxy will mainly contribute to devel- 
oping the massive galaxy's bulge component. The Hubble 
sequence then naturally appears as the result of the inter- 
play between cooling - which develops discs - and merging 
and disc gravitational instabilities - which develop bulges. 

Keeping track of the stellar content of each galaxy, as a 
function of age and metalicity, and knowing the galaxy's 
gas content and chemical composition, one can compute 
the (possibly extincted) spectral energy distribution (SED) 
of each galaxy. To thi s end, we use the Stardust model 
l|Devriendt et al.1 Il999l ) which predicts the SED of an ob- 
scured stellar population from the UV to the sub-mm. 



